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Magnetic model for A2CUP2O7 (A = Na, Li) revisited: 
one-dimensional versus two-dimensional behavior 
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J. Sichelschmidt,^ Yu. Skourski,'* G. Amthauer,^ and H. Rosner^^LO 

^ Department of Materials Engineering and Physics, Universitat Salzburg, 5020 Salzburg, Austria 

^ Max Planck Institute for Chemical Physics of Solids, 01187 Dresden, Germany 
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* Dresden High Magnetic Field Laboratory, Helmholtz-Zentrum Dresden-Rossendorf, 01314 Dresden, Germany 

We report magnetization measurements, full-potential band structure calculations, and micro- 
scopic modeling for the spin-1/2 Heisenberg magnets J42CUP2O7 (A = Na, Li) involving complex 
Cu-O-O-Cu superexchange pathways. Based on a quantitative evaluation of the leading exchange 
integrals and the subsequent quantum Monte-Carlo simulations, we propose a quasi-one-dimensional 
magnetic model for both compounds, in contrast to earlier studies that conjectured on a two- 
dimensional scenario. The one-dimensional nature of J42CUP2O7 is unambiguously verified by mag- 
netization isotherms measured in fields up to 50 T. The saturation fields of about 40 T for both Li 
and Na compounds are in excellent agreement with the intrachain exchange Ji ~ 27K extracted 
from the magnetic susceptibility data. The proposed magnetic structure entails spin chains with 
the dominating antiferromagnetic nearest-neighbor interaction Ji and two inequivalent, nonfrus- 
trated antiferromagnetic interchain couplings of about O.OlJi each. A possible long-range magnetic 
ordering is discussed in comparison with the available experimental information. 

PACS numbers: 75.30.Et, 75.50.Ee, 71.20.Ps, 75.10.Pq 



I. INTRODUCTION 

Quantum effects in magnetic systems have diverse 
implications in exotic ground statesi"— and interesting 
finite-temperature properties, such as spin transport!^ 
and multiferroic behavior j^i^ While transition-metal com- 
pounds offer nearly endless opportunities for finding vari- 
able spin lattices with different types of exchange cou- 
plings, complete understanding of the magnetic phenom- 
ena requires detailed and accurate information on the 
spin lattice and the energies of individual interactions. 
It is particularly important to establish the dimensional- 
ity of the system and the presence of magnetic frustra- 
tion. The latter may lead to degenerate ground states, 
while both are responsible for the strength of quantum 
fiuctuations. 

In complex crystal structures, the precise evaluation 
of the spin lattice remains a challenging problem. Sin- 
gle experimental observations on thermodynamic prop- 
erties or magnetic structure may lead to wrong conclu- 
sions regarding the nature of the spin lattice, as in the 
alternating-chain compound (VO)2P207 initially consid- 
ered as a spin ladder ji and the frustrated-spin-chain sys- 
tem LiCu202 that was also ascribed to spin laddersi^ 
Moreover, even the correct topology of the spin lattice 
does not guarantee the precise evaluation of the frustra- 
tion regime: see, for example, the studies of Li2VOX04 
with X = Si, GeJ^^— Generally, only a complete survey 
of the magnetic excitation spectrum or a comprehensive 
investigation of the ground state, thermodynamics, and 
electronic structure give unambiguous information on the 
underlying spin model. 

A particularly deceptive situation concerns the inter- 
pretation of the low magnetic ordering temperature T/v, 



as compared to the effective energy scale of exchange cou- 
plings, which is typically measured by the Curie- Weiss 
temperature 9. If T/v ^ |^|, long-range magnetic or- 
der is impeded by quantum fluctuations, but it is im- 
possible to decide a priori whether this effect originates 
from the low dimensionality, from the frustration, or 
from a combination of both. For example, a quasi-two- 
dimensional (2D) non-frustrated system should feature 
Tn/\&\ > 0.2,— while any lower value of T/v indicates 
sizable frustration. However, in a quasi-one-dimensional 
(ID) system Tn/\9\ may be as low as 0.01 even without 
frustration J^ A recent example of such an ambiguous sit- 
uation is (NO)Cu(N03)3 that was initially understood as 
a strongly frustrated quasi-2D system according to the 
magnetic susceptibility data and supposedly low T/vJ^ 
Subsequent electronic structure calculations did not find 
any signatures of the frustration, and rather showed the 
pronounced quasi-one-dimensionality that is a plausible 
reason for the low T^r.— While ab initio computational 
methods are an invaluable tool for studying complex ma- 
terials, independent experimental evidence is indispens- 
able. In the following, we show how both experimental 
and computational methods are separately used for eval- 
uating the dimensionality of the system. After the di- 
mensionality is established, simple criteria, such as the 
Tn/\6\ ratio, can be applied to analyze the frustration. 

We present our approach for the rather simple and 
non- frustrated, albeit controversial spin- 2 model com- 
pounds A2CUP2O7 {A = Li, Na). Both systems are 
antiferromagnetic (AFM) insulators and reveal appar- 
ent ID features of the crystal structure (Fig. H] and 
Sec. IIII[) . However, the magnetic dimensionality is hard 
to decide a priori, because both intrachain and inter- 
chain couplings involve complex Cu-0-O-Cu superex- 



change pathways that can not be assessed in terms of sim- 
ple empirical considerations, such as the Goodenough- 
Kanamori-Anderson rules.—"— The experimental data 
for Na2CuP207 reported by Nath et alt^ did not fit to 
the ID spin model, and showed better agreement with 
the 2D scenario. Later, Salunke et al^ performed elec- 
tronic structure calculations and argued that the leading 
couplings run along the structural chains, but the in- 
terchain couplings are stronger than in other quasi- ID 
compounds (e.g., in the chemically related K2CUP2O7 
presented in Ref. [21I ). and might be responsible for 2D 
features of the susceptibility at low temperatures. Here, 
we resolve the controversy and unequivocally establish 
the magnetic model of A2CUP2O7 by revising both ex- 
perimental and computational results. We confirm the 
quasi-lD scenario in an independent and more accurate 
electronic structure calculation, compared to the earlier 
work by Salunke et al^ We revisit the susceptibility data 
and present the results of high-field magnetization mea- 
surements that unambiguously evidence the quasi-lD na- 
ture of Na2CuP207 and of its Li analog Li2CuP207. 

The paper is organized as follows. In Section [III the ap- 
plied experimental and computational methods are pre- 
sented. Afterwards, in the third Section, the crystal 
structures are described and compared with similar com- 
pounds. In Section HVl the results of our theoretical in- 
vestigations and the experimental data are discussed. Fi- 
nally, the conclusions and a short outlook are given in 
Section El 



orbital code (FPLO) version 9.00-34^ For the scalar- 
relativistic calculations within the local density approxi- 
mation (LDA) , the Perdew-Wang parameterization of the 
exchange correlation potential^! was used together with 
a well converged k-mesh of 10 x 10 x 10 points. The strong 
electron correlations, only poorly described in LDA, were 
considered: (i) by mapping the magnetically active (par- 
tially filled) bands first onto a tight-binding (TB) model 



Htb = 






H.c. 



(1) 



and in a second step onto a single-band Hubbard model 
H = Htb + t^cff X)i"'it"'ii- ^"^ ^^^ strongly correlated 
limit, Ucs 3> tij, and for half- filling, well justified for 
undoped cupratcs, the low-energy sector of the Hubbard 
model may further be mapped onto a Hciscnberg model 



H 



7 , JijSiSj, 
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yielding the AFM contribution to the superexchange cou- 
pling constant J in second order as J^^^ = itfj/UcS- 
(ii) Alternatively, the effects of correlation were con- 
sidered in a mean-field way within the local spin den- 
sity approximation (LSDA)+C/ method {Ud = 7 ± 0.5 eV, 
Jd = 1 eV) J^i^i2£ For the double-counting correction, the 
around-mean-field approach^! was appliedi^ The mag- 
netic coupling constants Jy were obtained from energy 
differences between several collinear spin configurations. 



II. METHODS 

Bluish-colored powder samples of Na2CuP207 and 
Li2CuP207 were prepared through solid-state reaction 
technique by mixing NaH2P04-H20 (99.999% pure) or 
LiH2P04-H20 (99.99% pure) and CuO (99.99% pure) in 
appropriate molar ratios. The stoichiometric mixtures 
were fired at 800 °C and 750 °C for 60 hours each, with 
one intermediate grinding and pelletization to achieve 
single-phase Na2CuP207 and Li2CuP207 compounds, 
respectively. The purity of the samples was confirmed 
by x-ray diffraction (Huber G670 camera, CuKqi radia- 
tion, ImagePlate detector, 29 = 3 — 100° angle range). 

Magnetic susceptibility was measured with the com- 
mercial MPMS SQUID magnetometer in the temperature 
range 2 — 380 K in applied fields up to 5 T. High-field data 
were collected at a constant temperature of 1.5 K using 
pulsed magnet installed at the Dresden High Magnetic 
Field Laboratory. Details of the experimental procedure 
can be found elsewherei^ 

Electron spin resonance (ESR) spectra were measured 
at room temperature using the X-band frequency of 
9.4 GHz. The spectra were fitted as powder average of 
several Lorentzian lines reflecting different components 
of the anisotropic g-tensor. 

The electronic structure was calculated using the den- 
sity functional theory (DFT)-based full-potential local- 



Quantum Monte-Carlo (QMC) simulations were per- 
formed using the looper— and dirloop_sse (directed 
loop in the stochastic series expansion representation)^ 
algorithms of the software package ALPSi^ The magne- 
tization was simulated for the N = 400-site chain with 
periodic boundary conditions. To evaluate the ordering 
temperature T/v , we considered the full three-dimensional 
(3D) spin lattice entailing interchain couplings, and 
calculated the Binder ratio^ B = {Mf)/{M^)'^ for 
the staggered magnetization Mg on the cluster size N. 
We performed a series of simulations starting with a 
N = 640-sites cluster and consequently increasing it up 
to Af = 40960 sites. In the temperature range 0.1 < 
T/Ji < 18.0, we used periodic boundary conditions, 
20 000 sweeps for thermalization and 200 000 sweeps af- 
ter thcrmalization. The resulting statistical errors (below 
0.1 %) arc negligible compared to the experimental error 
bars. Magnetic susceptibility of a Heisenberg model on 
a square lattice was computed on the iV = 20 x 20 finite 
lattice, using 30 000 sweeps for thermalization and 300 
000 sweeps after thermalization. 



III. CRYSTAL STRUCTURE 

Na2CuP207 crystallizes in two dissimilar mono- 
clinic structurcSf2S In the following, we focus on 
the /3-polymorph that is isostructural to Li2CuP2 07. 
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FIG. 1. (Color online) Crystal structure of Na2CuP207 and 
the leading magnetic coupling pathways. The indices ic 
and ip denote interchain and interplane transfer, respec- 
tively. Apart from substitutions of Na by Li, the structure 
of Li2CuP2 07 shows only minor differences. 



/3-Na2CuP207, further referred as Na2CuP207, has the 
space group C2/c with the lattice parameters a = 
14.728 A, b = 5.698 A, c = 8.067 A, and P = 115.15°^ 
The crystal structure entails chains stretched along the c 
direction. Each chain consists of CUO4 plaquettcs linked 
via two corner-sharing PO4 tetrahedra. Neighboring pla- 
quettcs arc tilted toward each other by an angle of about 
70° (Fig. [1]). Compared to a planar arrangement of the 
plaquettes in other quasi-lD Cu"*"^ phosphates, such as 
K2CUP2O7 exhibiting a pronounced quasi-lD magnetic 
behavior)^ the tilting could enhance the interchain cou- 
plings and might be responsible for the proclaimed quasi- 
2D magnctismJ^i^S On the other hand, CUSC2O5 with a 




-4 -2 

Energy (eV) 

FIG. 2. (Color online) LDA electronic density of states 
and partial density of states at the Fermi level (inset) 
for Na2CuP207 arising from the two well separated bands 
(see Fig. [3|. Oljpa) denotes the contribution from the a- 
antibonding 2p-orbitals of the O atoms being part of the pla- 
quettes. O(term) denotes the contributions from O-atoms be- 
ing not part of the plaquettes. The results for Li2CuP207 are 
very similar (see Fig. [3]). 



2.5 
2 

1.5 

1 

0.5 



^ 30 
I 25 
S 20 

% 10 
w 5 






-0.2 -0.1 0.1 0.2 

Energy (eV) 




FIG. 3. (Color online) Comparison between the LDA band 
structures of Na2CuP207 and Li2CuP207. The inset shows 
the DOS of both compounds at the Fermi level. 



similar tilting angle as in Na2CuP207 was clearly shown 
to be of ID type,— although the distances between the 
6c-planes are significantly smaller and even the chains 
within the 6c-plane are somewhat closer together than 
in Na2CuP207. Accordingly, it could be expected from 
structural considerations that Na2CuP207 is a quasi-lD 
system. 

Li2CuP207 is isostructural to /3-Na2CuP207 with the 
unit cell parameters a = 15.336 A, b = 4.8733 A, c = 
8.6259 A, and P = 114.79°.-^>^ The tilting angle is about 
90°, leading to shorter distances between the chains 
within the 5c-plane, where in particular the ii^ic path is 
significantly shorter. The distances between the planes 
arc only slightly smaller than in Na2CuP207. Since in 
these compounds the plaquettes are linked via the phos- 



TABLE I. Transfer integrals (ti) and AFM contributions to 
the excliange integrals (J^ ^^ = At\/Uc«) for Na2CuP207 and 
Li2CuP207, as obtained from full-potential LDA calcula- 
tions. The hoppings by Salunke et al^ were calculated for 
Na2CuP2 07 within the ASA (see text for details). 
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phate tetrahcdra, the empirical Goodcnough-Kanamori- 
Anderson rule&i^"— cannot be applied for describing the 
NN superexchange as a function of the tilting angle, and 
an elaborate analysis of the electronic structure is re- 
quired. 



IV. RESULTS AND DISCUSSION 
A. DFT calculations 

We start with the computational analysis of 
A2CUP2O7. The DFT calculations of the band struc- 
ture and the electronic density of states (DOS) within 
the LDA yield a valence band width of about 9 eV and 
9.5 eV (Fig.© for NasCuPzOy and LizCuPzOr, respec- 
tively, typical for cuprates. The band structures of the 
two compounds are very similar and show well separated 
bands crossing the Fermi level (Fig. [3|) somewhat nar- 
rower in the case of Li2CuP207. Typical for cuprates, 
these bands are formed by cr-antibonding linear combi- 
nations of Cvi{3dx2_y2) and 0{2p) orbitals (Fig. [21 inset). 
The orbitals are denoted with respect to a local coordi- 
nate system, where for each plaquette one of the Cu-0 
bonds and the direction perpendicular to the plaquette 
are chosen as x- and z-axes, respectively. 

Since the two bands at the Fermi level corresponding 
to the two Cu atoms per unit cell are well separated from 
the rest of the valence band, the magnetic properties can 
be described using an effective one-orbital model. Al- 
though LDA yields a wrong metallic ground state due 
to the poor description of strong electronic correlations, 
the LDA bands can be used as an input for an effective 
Hubbard model that introduces the missing correlations 
and yields the AFM part of the exchange, as described 
in Sec. [ill A crucial parameter in this computational 
scheme is the LDA bandwidth W determined by individ- 
ual hoppings ti. The bandwidth may be affected by the 
approximations introduced in the DFT calculation. Par- 
ticularly, the atomic spheres approximation (ASA) for 
the potential is known to underestimate the bandwidth 



in Cu"*"^ compounds 41 Since the previous computational 
work by Salunke et al^ was based on the ASA calcula- 
tion, the insufficient accuracy of the band structure could 
lead to a sizable error in the hoppings and exchange cou- 
plings, thus hindering the reliable evaluation of the di- 
mensionality. 

To avoid the problems of the ASA, we performed full- 
potential LDA calculations that produce a highly ac- 
curate estimate of the LDA bandwidth and individual 
hoppings 1^ The hoppings are obtained from Wannier 
functions with the Cu 3dx2-y2 orbital character^ and 
perfectly reproduce the LDA bands, as demonstrated in 
Fig. [5] In Table U we compare our results for both 
Li2CuP207 and Na2CuP207 with the previous ASA data 
by Salunke et al^ for Na2CuP207. In the ASA calcula- 
tions, the bandwidth W (Fig. |4]) and the intrachain hop- 
ping ii ~ W/A are indeed underestimated for about 25 %, 
which is similar to our previous comparative study of the 
ASA and full-potential calculations for Sr2Cu(P04)2'^ 

Applying the effective on-site Coulomb repulsion 
Ucs ~ 4.5 oV j2i,28.39 ^g evaluate AFM contributions to 
individual exchange couplings using the second-order ex- 
pression J^F^ = 4if /t/cff. (Table|Il)4a Since the spin lat- 
tice does not contain triangular units;^ the higher-order 
corrections^^iM are proportional to tj/U^ff and rather 
small compared to both intrachain and interchain cou- 
plings Jj^^^. For example, the fourth-order correction 
to the bilinear exchange is — 24tf/C/cfr — —0.1 K and the 
leading four-spin term is SOtft'^/U^g ~ 0.03 K, where 
ta = t2^ic,tipi^ Both fourth-order terms are comparable 
to or even smaller than typical dipolar interactions in a 
spin-i compound.— 

While individual hoppings derived from the full- 
potential and ASA calculations are rather different (note, 
e.g., the different signs of ^2,10), the resulting microscopic 
scenario is very similar. The interchain couplings J2,ic 
and Jip are about two orders of magnitude smaller than 
the intrachain coupling Ji. Although Salunke et al^ 
claim that "the interchain interactions in Na2CuP207 
can not be neglected" , our estimates of J^™ and J^^^ 
suggest a pronounced ID character of ^2CuP207. In 
fact, Na2CuP207 features comparable couplings along 
the two interchain directions (J^™ — J:^^^), and 
should be considered either as a quais-lD or as a three- 
dimensional (3D) system, in contrast to the 2D behavior 
proposed in the experimental study.— 

While the discrepancy between the computational and 
experimental results is actually resolved on the experi- 
mental side (Sec lIVB]) . we also check for possible short- 
comings of our calculations, and supplement our model 
estimate of J^^^ with the full exchange couplings Jt ob- 
tained from LSDA-hC/ calculations. Such full exchange 
couplings were not evaluated by Salunke et al^ The 
FM contributions considered in the LSDA-I-C/ calcula- 
tions reduce the NN coupling constant Ji to 35 ± 4K 
and 34 ± 4K for Na2CuP207 and Li2CuP207, respec- 
tively (compare to J^^^ in Table |T|. Interchain and 
intcrplane couplings were estimated to be smaller than 
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FIG. 4. (Color online) Comparison between the density of 
states at the Fermi level calculated using a full-potential code 
(this study) and an ASA code (Ref.l20[l. 




FIG. 5. (Color online) The LDA bands at the Fermi level and 
the band dispersions calculated with a tight-binding model 
for Na2CuP207 (upper panel) and Li2CuP207 (lower panel). 



0.2 K. Thus, the LSDA-hC/ results further support the 
ID scenario. Sunilar ratios of the interchain (Jic/Ji) 
and the interplane (Jip/Ji) couplings were found in sev- 
eral other quasi-lD systcmS i'^^''^^ In conclusion, unlike in 
previous studics22. we find no indications for the 2D be- 
havior neither for Na2CuP207 nor for Li2CuP207 based 
on electronic structure calculations. 



B. Experimental data 

After establishing the quasi-lD microscopic scenario, 
we re-consider the experimental data to check whether 
the previous conclusion on the 2D magnetic behavior- 
was well-justified. The conjecture about the 2D mag- 
netism made by Nath et alt^ was essentially based on the 
better agreement with the experimental data (nuclear- 



magnctic-rcsonance shift as a measure of intrinsic mag- 
netic susceptibility), in comparison to a ID model. How- 
ever, our precise QMC simulations for both ID and 2D 
models challenge this interpretation. 

Magnetic susceptibility (x) of the Li and Na com- 
pounds is very similar and shows the maximum around 
16 K, as typical for low-dimensional magnets (Fig. [5]). 
The increase in x below 7 K is likely related to the Curie- 
like contribution of paramagnetic impurities. We fit the 
data assuming simplest ID and 2D coupling topologies: 
a NN chain and a square lattice. Reduced magnetic sus- 
ceptibility X* was simulated using QMC for large finite 
lattices of iV = 400 and A^ = 20 x 20 spins, for the chain 
and the square lattice model, respectively. The simulated 
X*{T /[kB Ji\) dependency was fitted to the experimental 
curves using the expression: 



xiT) 
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(3) 



where the fitted parameters are the exchange cou- 
pling Jl, the Lande factor g, the extrinsic (impurity 
and defect) paramagnetic contribution Cimp/T, and the 
temperature-independent term xo that accounts for the 
core diamagnetism and van Vleck paramagnetism. 

The resulting fits are shown in Fig. [SJ with the fit- 
ting parameters listed in Table HIl For both Na2CuP207 
and Li2CuP207, the square-lattice model apparently 
fails to describe x{T) below 10 K. In contrast, the NN 
chain model provides a much better description of the 
data down to 2K. This contradicts the previous conclu- 
sion by Nath et al.^ although our susceptibility data 
for Na2CuP207 are very similar to those previously re- 
ported. The discrepancy originates from different model 
expressions used in our analysis. Nath et alJ^ utilize 
the high-temperature expressions that are valid down to 
T ~ Ji/2 ~ 14 K only. Above 14 K, the ID and 2D fits 
nearly coincide, hence neither model can be chosen un- 
ambiguously. The QMC simulations enable the precise 
evaluation of the susceptibility down to the lowest tem- 
perature of 2 K accessed in the present experiment. The 
low-tcmpcrature susceptibility data clearly favor the ID 
model, while rule out the 2D description of J42CUP2O7. 
Note that the experimental estimates of Ji ~ 27K within 
the ID model are in reasonable agreement with Ji ~ 35 K 
computed by LSBA+U. 

Four adjustable parameters of the susceptibility fit 
might lead to an ambiguous solution. Therefore, we 
verified the fitted g-values by an ESR measurement. 
The room-temperature ESR spectrum of Li2CuP207 was 
well described by two components of the g-tcnsor. ^n = 
2.064(2) andg_L = 2.198(2). In caseof Na2CuP207, three 
components were required: g^ = 2.217(5), g^ = 2.152(3), 
and g^ = 2.085(6). Both sets of the g- values result in 
the powder-averaged g ~ 2.15(1), in good agreement 
with the fitted values of 2.12 and 2.19 for Li2CuP207 
and Na2CuP207, respectively (see the fits for the chain 
model in Table [TTI'^ The three g- values resolvable in 
the ESR spectrum of Na2CuP207 might be related to 



10- 



'o 5- 



10 




r(K) 



FIG. 6. (Color online) Fits of the magnetic susceptibil- 
ity of Na2CuP2 07 (top) and Li2CuP2 07 (bottom). Mag- 
netic susceptibility was simulated using QMC for a Heisen- 
berg nearest-neighbor one-dimensional (chain) and two- 
dimensional (square lattice) models. Insets: simulated mag- 
netization curves on top of the experimental high-field mag- 
netization data. For the scaling, the parameters from the fits 
to x{T) were used (see text for details).— 



the distorted CUO4 plaquette (Cu-0 bond lengths of 
1.91 — 1.96 A) compared to the nearly regular CUO4 pla- 
quette in LiaCuPzOr (four Cu-0 bonds of 1.93 A). 

Although the susceptibility fits are already a good ev- 
idence for the quasi-lD magnetism, one might still argue 
that the ID model is solely chosen based on the low- 
temperature data that contain sizable impurity contri- 
butions. Moreover, the ESR results do not allow to dis- 
criminate between the marginally different (^-values ob- 
tained in the chain and square- lattice fits (Table HI)) . To 
underpin our conclusions, we measured magnetization 
isotherms up to 50 T, and observed the magnetic sat- 
uration of the AaCuPaOy compounds. The saturation 
field (Hs) is an independent measure of the exchange 
couplings. The comparison between Hs from the magne- 
tization curve and Ji from the susceptibility fit is a sim- 
ple and efficient test for the validity of a spin model4i 
Experimental data (insets to Fig. |6]) show that the mag- 
netization of both Li and Na compounds saturates above 
40 T. The upward curvature below Hg is a typical fea- 
ture of low-dimensional spin systems, and is generally 



TABLE II. Fitting parameters for the ID and 2D descrip- 
tion of the magnetic susceptibility of A2CUP2O7: the intra- 
chain exchange coupling Ji, the g- value, the temperature- 
independent contribution xo, and the Curie-like paramagnetic 
impurity contribution Cimp. 
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ascribed to quantum effects.— 

To enable the proper comparison between the exper- 
imental and simulated magnetization, we subtract the 
paramagnetic impurity contribution according to Cimp 
obtained from the susceptibility fit (2 — 5% of the para- 
magnetic impurity depending on the model and com- 
pound). In Fig. ini we juxtapose the experimental data 
and the simulations^ for both ID and 2D models using 
the parameters from Table HH as derived from the suscep- 
tibility measurements. The excellent agreement between 
the ID scenario and the magnetization data verifies our 
suggestion on the quasi-lD model, while the 2D model 
apparently fails to reproduce the high-field data. Note 
that our conclusion is robust with respect to the impurity 
contribution, because the deviations from the 2D model 
are observed at high fields, where the paramagnetic con- 
tribution is saturatued and field-independent. 

We have shown that the high-field magnetization data 
are an efficient tool for distinguishing between different 
spin models and even determining the dimensionality of 
the system. It is instructive to consider how the magne- 
tization curve resolves the ambiguity of the susceptibility 
fit. In simple spin models, the position of the susceptibil- 
ity maximum uniquely determines the exchange coupling 
Ji. For example, Tmax/^i — 0.64 and Tmax/^i — 0.95 for 
the uniform chain^ and square lattice, respectively. This 
explains the 30% difference in the respective estimate of 
Ji (Table |Tl|. The saturation field is directly related 
to the energies of different magnetic states, according to 
Hs = 2JiA:b/(5/^b) and Hs = 4Ji/cb/((7//_b) for the ID 
and 2D cases, respectively. Therefore, the same satura- 
tion field Hs leads to the 50% difference in Ji . 

The dissimilar dependences of T^ax and Hs on Ji pro- 
vide a key to discriminate between the ID and 2D mag- 
netism using the combination of susceptibility and high- 
field magnetization data. The sizable difference between 
the saturation fields of the ID and 2D models ensures 
that this method can be applied even to weak couplings 
on the order of 5 K that produce about 1 T difference in 
Hs- Systems with larger couplings are, of course, easier 
to evaluate, although the couplings above 50 K shift the 
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FIG. 7. (Color online) QMC simulations of the Binder ra- 
tio B{T) = {Mf}/{M^f of staggered magnetization M^ for 
Jic/Ji =0.05 (left panel) and Jic/Ji =0.01 (right panel). The 
finite lattice dimensions are given in the legend. The ordering 
temperatures Tjv/Ji are highlighted. 



saturation field above the limits of present-day facilities. 
Our simulations shown in the insets to Fig. |6] suggest that 
even the part of the curve right below Hg (at 30 — 35 T 
in the present case) may be sufficient to discriminate be- 
tween the ID and 2D models. However, the data at lower 
fields do not contain the necessary information, because 
below 30 T the magnetization curves for the ID and 2D 
systems nearly coincide. 



C. Long-range ordering 



In Sec. lIVBl we have shown that thermodynamic prop- 
erties of A2CUP2O7 above 2 K are well reproduced by the 
purely ID model of the uniform spin chain. The actual 
system is, however, only quasi- ID, because any chem- 
ical compound necessarily entails both intrachain and 
interchain couplings. In A2CUP2O7, the non-frustrated 
interchain couplings Jip and Ji^ic (Li2CuP207) or J2,ic 
(Na2CuP2 07) should manifest themselves at low temper- 
atures and eventually lead to the long-range magnetic 
ordering. No definitive signatures of the magnetic order- 
ing were reported by Nath et al.,— according to heat- 
capacity and nuclear magnetic resonance (NMR) mea- 
surements down to 2 K, although the broadening of the 
NMR line could indicate the approach to the magnetic 
transition below 4K. Our susceptibility measurements 
down to 2 K did not reveal any anomalies in both Li 
and Na compounds cither, thus suggesting the lack of 
the long-range order down to 2K. 

The fact that the A2CUP2O7 compounds do not or- 
der down to 2 K {T/Ji ~ 0.08) is an additional argument 
against the quasi-2D scenario, because the non-frustrated 
quasi-2D system orders at T/Ji > 0.2.ii^ The low Tn is, 
however, in excellent agreement with the quasi- ID sce- 
nario and indirectly supports our conclusions. To eval- 
uate T^v, we considered the realistic 3D magnetic model 



of J42CUP2O7 with the intrachain coupling Ji as well as 
the interchain couplings J2,ic and Jip represented by a 
single effective interchain coupling Ji^. We utilize the 
scaling property of the Binder ratio (B) to be indepen- 
dent of the cluster size at T/v^ Therefore, T^r is deter- 
mined as the crossing point of B{T) curves calculated 
for clusters of different size. The typical scaling is shown 
in Fig. [71 where we consider the cases of Jic/Ji — 0.05 
and Jic/Ji = 0.01. For the former case, T/v/Ji ~ 0.13 
equals to 3.5 K and should be visible in the present ex- 
periments. Taking the weaker interchain coupling as the 
realistic estimate for Na2CuP207 (Table |T|, we arrive at 
T/v/Ji ~ 0.032, which is as low as 0.9 K, well below the 
temperature range studied in our work. In Li2CuP207, 
a similar Tn is expected, because Jiac is comparable to 
J2.ic in the Na compound. 

The anticipated Tn/Ji ^ 0.032 in Na2CuP207 is not 
the lowest transition temperature reported for quasi-lD 
magnetic systems. For example, Sr2Cu(P04)2 lacks the 
long- rang e magnetic order down to T/v/ J = 4.5 x lO"'^ 
(Ref. [50 and Isil ). This anomalously low T/y is, how- 
ever, facilitated by the frustrated nature of interchain 
couplingsi^ In other spin-chain compoundsfi^ the in- 
terchain couplings are often anisotropic, i.e., the inter- 
chain couplings in the plane are stronger than along the 
third direction. In contrast to all these perplexing scenar- 
ios, Na2CuP207 reveals an unusually simple geometry of 
the interchain couplings with similar interactions along 
a and b. This coupling regime conforms to the quasi- 
lD mode l typ icall y co nsidered in theoretical studies (e.g., 
Refs. I13I. [53. andlSSh. Therefore, it may be interesting 
to explore the magnetic transition in Na2CuP207 exper- 
imentally. 

Experimental evaluation of T/v and long-range mag- 
netic order will, on one hand, verify theoretical results 
for weakly coupled spin chains, and on the other hand 
challenge our microscopic magnetic model. One efficient 
test is the magnetic structure that should feature an- 
tiparallel spins along a and c, yet parallel spins along b 
owing to the "diagonal" interchain coupling J2,ic (Fig.[T]). 
The anticipated propagation vector is k = 0, because the 
crystallographic unit cell contains two Cu sites along a 
and c, while the order along b is FM. 



V. SUMMARY AND OUTLOOK 

Motivated by suggestions of a quasi-2D behavior of 
Na2CuP207; exhibiting a crystal structure similar to well 
known quasi- ID-compounds, we have reinvestigated its 
microscopic model along with the previously unexplored 
Li2CuP207 compound. To this end, magnetization mea- 
surements, full-potential DFT calculations, and quantum 
Monte Carlo simulations have been applied. We inde- 
pendently confirmed the anticipated quasi-lD scenario 
by the experimental data and by numerical estimates of 
individual exchange couplings. The computed intrachain 
couplings arc in excellent agreement with the experiment. 



The previous conjecture on the quasi-2D magnetic behav- 
ior is ascribed to the incomplete analysis of the magneti- 
zation data. 

The microscopic model of Na2CuP207 is a rare exam- 
ple of a "regular" quasi- ID system with similar couplings 
along the two interchain directions. Therefore, this com- 
pound is an excellent prototype material for the simplest 
model of weakly coupled spin chains that was widely 
studied theoretically. A further experimental study of 
the anticipated long-range ordering at Tn ~ 0.9K may 
be insightful, as outlined in Sec. IIVCI 

Another interesting aspect of the ^2CuP207 com- 
pounds are magnetostructural correlations for the su- 
perexchangc in Cu"''^ phosphates and other transition- 
metal compounds containing polyanions. The intrachain 
couplings in Li2CuP207 and Na2CuP207 {J\ ~ 27 K) 
are much smaller than in K2CUP2O7 ( Ji ~ 141 K)2i and 
Sr2Cu(P04)2 (Ji ^ 150 K)^ with similar chains of CUO4 
plaquettes bridged by PO4 tetrahedra. The difference 
could be ascribed to the alteration of the chain geometry: 
both K2CUP2O7 and Sr2Cu(P04)2 feature the coplanar 
plaquettes, while in A2CUP2O7 the neighboring plaque- 
ttes form an angle of 70 - 90° (Sec. [ml). The buckhng 
impedes the overlap of the pa oxygen orbitals on the O- 



O edges of the PO4 tetrahedra, thereby reducing the ex- 
change. However, this simple explanation does not hold 
for CuSe205, where the buckling angle of 64° resembles 
that of y42CuP207, but the strong intrachain exchange 
of J ~ 157 K rather reminds of the planar geometry in 
K2CuP207.'^ A plausible explanation is the different role 
of the non-magnetic P^^ and Se"*"^ cations, and a further 
study of this matter for a broader compound family could 
be insightful. 
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